Vortex molecules in coherently coupled two-component Bose-Einstein condensates 
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A vortex molecule is predicted in rotating two-component Bose-Einstein condensates whose in- 
ternal hyperfine states are coupled coherently by an external field. A vortex in one component 
and that in the other are connected by a domain wall of the relative phase, constituting a "vortex 
molecule" , which features a nonaxisymmetric (pseudo)spin texture with a pair of merons. The bind- 
ing mechanism of the vortex molecule is discussed based on a generalized nonlinear sigma model and 
a variational ansatz. The anisotropy of vortex molecules is caused by the difference in the scattering 
lengths, yielding a distorted vortex-molecule lattice in fast rotating condensates. 
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Topological defects appear in cross-disciplinary sub- 
fields of physics as long-lived excitations, constrained by 
the topology of the order parameter [ij. A prime ex- 
ample is quantized vortices which play a key role in un- 
derstanding of superfluidity 2]. When a system has a 
multicomponent order parameter, it is possible to excite 
various exotic topological defects which have no analogue 
in systems with a single-component order parameter. 

An atomic-gas Bose-Einstein condensate (BEC) offers 
an ideal testing ground to investigate such topological 
defects, because almost all parameters of the system 
can be controlled to the extent that state engineering 
is possible. Because alkali atoms have a spin degree of 
freedom, multicomponent BECs can be realized if more 
than one hyperfine spin state is simultaneously populated 
Pl yi ' A quantized vortex and a vortex lattice in BECs 
have been created experimentally by several techniques 
Himii' The methods reported in Ref. [IIH utilized 
internal degrees of freedom of BECs, creating unconven- 
tional vortices described by multicomponent order pa- 
rameters. The structure of single vortex states in systems 
with multicomponent order parameters was investigated 
in Refs. [EllflllJ- In addition, it has been predicted 
theoretically that fast rotating two-component BECs ex- 
hibit a rich variety of unconventional vortex structures 

[Ulillii. 

In this Letter, we study the vortex structure of ro- 
tating two-component BECs whose internal states are 
coupled coherently by an external driving field. This cou- 
plin g ca n be achieved experimentally as reported in Refs. 
[ij, ll5| , where Rabi oscillations between the two com- 
ponents were observed. If the strength of the coupling 
drive is increased gradually from zero and its frequency 
is gradually ramped to resonance, one can obtain a sta- 
tionary state with a nearly equal-weight superposition of 
the two states ^^. Here, we study stationary states of 
two-component BECs with an external rotation as well 
as internal coupling. Combination of these two effects en- 
ables us to explore a new regime of rich vortex structures 
beyond the conventional binary system jSlU^lldilld) ^^^ 



two components interact not only through their mean- 
field interactions but also through the relative phase of 
the order parameters. We find that such two-component 
BECs exhibit unique vortex structures; a vortex in one 
component and that in the other form a stable vortex- 
antivortex pair |l7| . i.e., a "vortex molecule" , bound by a 
domain wall in the space of relative phase. The psudospin 
representation of the molecule reveals a nonaxisymmet- 
ric spin texture with a pair of merons |l8l lltA |20j . The 
internal coupling controls the binding force of the vortex 
molecule and affects significantly the lattice structure of 
vortices in fast rotating two-component BECs. 

We consider two-component BECs of atoms with mass 
m residing in two hyperfine states in a rotating frame 
with rotation frequency ft = ilz. The BECs are de- 
scribed by the condensate wave functions ^i {i = 1,2). 
The equilibrium state of the system is obtained by min- 
imizing the energy functional 

i?[*i,*2]=yrfr|^ U*h,^, + ^\^,\ 



In this work, we confine ourselves to the two-dimensional 
problem. Here, we measure the length and energy in 
units of 6ho = \/h/mu} and hcu, respectively, with a fre- 
quency w of a harmonic trap; thus hi = — V^/2 -I- r^/2 — 
riLz. The atom-atom interactions are characterized by 
Ui for intracomponent and M12 for intercomponent inter- 
actions. The last term in Eq. |^ describes a coher- 
ent coupling induced by an external field, which allows 
atoms to change their internal state coherently [ij, ll5j . 
Here wr is the Rabi frequency and A is the detuning 
between the external field and the atomic transition; in 
the following we set A = for simplicity. In this case, 
the total particle number (per unit length along the z 
axis) N = /dr{|vl/ip-|-|vl/2p} = J drpT is conserved. By 
renormalizing the wave function as V&j -^ \/iV*I'i /&ho, the 
atom-atom interactions are written as Ui = AiraiN and 
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FIG. 2: (a) The cross sections of the condensate density along 
the y = hne (|^ip: sohd-curve, |^2p: dashed-curve, px: 
dotted-curve). (b) The total energy as a function of A (the 
size of the meron pair) for f2 = 0. The inset shows the sep- 
aration 2dm between two vortices as a function of the Rabi 
frequency wr for f? = 0.15. The solid curve represents the 
result obtained by the variational analysis with Eqs. and 
@ (dm = 2 A) and the dots the numerical result. 



FIG. 1: (a) The profile of the density |*ip (color-scale plot) 
and j*2p (line contours) with fi = 0.15 and tjR = 0.02. (b) 
The gray-scale plot of the relative phase <^(r) — 9i — 62. Ar- 
rows show the direction of the circulation in the space of rel- 
ative phase around the vortices, (c) The vectorial represen- 
tation of the pseudospin S — x^x/S projected onto the x — y 
plane. The locations of the defects are indicated by (with 
Sz = +z/2) and (with 5*^ = -z/2). (d) The contour plot 
of the topological charge density g(r) (the largest value at the 
center) and the vectorial plot of the VcH-field. 



M12 — 4:T:ai2N with the corresponding s-wave scattering 
lengths Ui and ai2- 

The crucial effect of the internal coherent coupling 
is to induce an effective attractive interaction between 
the two components 16]. Putting ^^ — |^i|e'^', 
we can rewrite the last two terms in Eq. ([[J as 
-2cjR,|vl/i(r)||^I'2(i")| cos0(r), where 0(r) = 6*1-02 is the 
relative phase. To decrease it, must decrease and the 
amplitudes of the two components must increase in the 
overlapping region. 

We first assume ui = U2 = U12 — 1000 (a more general 
case will be discussed later) to discuss the vortex state 
with the external coupUng and an appropriate rotating 
drive to ensure stabilization of only one vortex in each 
component. The minimization of the energy functional 
(Q is done numerically by using the conjugate gradient 
method. Figure ^a) shows the equilibrium solution for 
wr = 0.02 and fl — 0.15. Each component has one off- 
centered vortex which shifts from the other to reduce the 
overlapping area [see Fig. Efa)]. For finite wr. this non- 
axisymmetric vortex state is always energetically lower 
than the axisymmetric one which was observed in Ref. 
[5| , in which one circulating component surrounds a non- 
rotating core of the other '^. Also, the optimized relative 



phase shows an unique structure as shown in Fig. Q^b). 
Here, the central region is characterized by the config- 
uration of a vortex-antivortex pair; the vortex cores are 
connected by a branch cut of the relative phase with a 27r 
phase difference |21|. Thus, the two vortices attract each 
other, forming a bound pair, i.e., a "vortex molecule". 
As wr increases, the size of the pair decreases as seen 
in Fig. Ela). Beyond ojr ~ 3.0 the separation vanishes, 
where the locations of the density nodes overlap despite 
the intercomponent repulsive interaction. 

An insight into the vortex molecule can be gained when 
we describe the two-component BECs in terms of a pseu- 
dospin |lll Il5| . By introducing a normalized spinor 
X = [Xi(i"),X2(r)]'^ with |xip + |x2p = 1 and writ- 
ing '^i = v^pT(r)Xi(r), we can define the spin vector 
as S(r) = xcrxl^ with the Pauli matrix <t; the modulus 
of the total spin is |S| = 1/2. The vectorial representa- 
tion of the spin density (projected onto the x — y plane) is 
shown in Fig. ^c). The spins are oriented in the x direc- 
tion everywhere except in the central domain-wall region 
where they tumble rapidly by 27r. There exist two points 
corresponding to the locations of vortices at which S is 
parallel to the z-axis. The spin field around the singular- 
ity with S — +z/2 (S — —z/2) has a radial (hyperbolic) 
distribution, having (Sx,Sy) oc (— x, — y) (oc (x, — y)). 
This texture is known as a "radial-hyperbolic" pair of 
merons l22l | , which has been discussed not only in the 
study of topological defects in superfiuid ^He jlj| and a 
double-layer quantum Hall system 19] but also in the 
semiclassical model of color confinement in QCD |20|. 

To discuss the properties of a vortex molecule quanti- 
tatively, we rewrite Eq. ^ in terms of the pseudospin 



variables as 
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where we have introduced the effective velocity field in- 
duced by spin textures 

^® ^' :{SyVS^-s^\iSy) (3) 



Voff 



(5j. + S ) 



with the total phase 8 = 6'i + 6*2, and three coupling 
constants Cq = {ui + U2 + 2ui2)/4, ci = wi — U2 and 
C2 = Ui + U2 ~ 2wi2. Equation ^ describes a gener- 
alized nonlinear sigma model for two-component BECs, 
which allows us to obtain the solutions for topologi- 
cal excitations such as "skyrmions" or "merons" |l9f . 
They are characterized by the topological charge density 
g(r) = e^^S • (9^S) x (c'^S)/7r = (V x Voff)/27r, whose 
spatial integral is the topological charge. Figure ^d) 
shows the profile of q(r) together with the vector plot of 
the Veff -field. Contrary to a conventional vortex, |veff| 
vanishes at the center, refiecting a coreless vortex with- 
out a density dip in the total density. The topological 
charge of this solution is exactly unit. The axisymme- 
try of both profiles is reflected to the SU(2)-symmetry 
associated with ci — C2 = 0. Then, the structure of the 
meron-pair is equivalent to an axisymmetric vortex state, 
i.e., a skyrmion jllj . where both states are connected via 
overall spin rotation with 90 degrees [22|; according to 
our defenition of S, the spin of the skyrmion points down 
along the z-axis at the center and rolls up continuously 
as it goes outward if the ^i component has a vortex. The 
internal coupling works as a transverse "magnetic" field 
that aligns the spin along the a;-axis. Therefore, turning 
on iOB. makes a skyrmion split into two merons with the 
spin texture shown in Fig. ^c) 19|. 

The pseudospin picture gives us an physical interpreta- 
tion of vortex-molecule binding by means of the internal 
coupling. When ci = C2 ~ 0, the pseudospin profiles of 
the meron-pair can be represented by the form 
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where A and a are the 



variational parameters that characterize the size of the 



the molecular size d?,. 



meron-pan, i.e., tiie molecular size dj^e""™ — 4A (see 
the top of Fig. |21 (a) for the definition of dm)- The 
problem becomes simpler by assuming a = so that 
Eq. I0J reduces to a solution of the classical nonlin- 
ear sigma model ll9|. Substituting S into Eq. (0) 
and putting = tan^^ J^„. + tan~^ _^„. , we obtain 
E = \ dr[(VVp^) 2/2 -f FeffPT + cop|/2] with 



Kff = Tr 
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FIG. 3: The contour plots of the topological charge density 
g(r) and the vectorial plots of the Veff-field for f2 = 0.15, ojr = 
0.02, Co = 1000 and (a) ci = 20 (antiferromagnetic), and (b) 
C2 = —20 (ferromagnetic). The inset shows the corresponding 
density profiles. 



To see the (jjR,-dependence on the molecular binding, we 
calculate the total energy £ as a function of A in the case 
without rotation |2J|; the result is shown in Fig. El^b). 
Here, for a given A(= dm/2) we calculate numerically 
Pt that minimizes E. For wr = the energy decreases 
monotonically with A, which implies the repulsive inter- 
action between two merons. Its dominant contribution 
comes from the second and third term of Eq. |(2Jl, which 
are the gradient energy of the pseudospin and the kinetic 
energy of the Veff-field. Since the total density /Ot is fixed 
as the Thomas-Fermi profile with an inverted parabola 
for the moderate values of A except A < 0.15, where the 
vortex core appears in pT as seen in Fig. I2a), the energy 
contribution of the other terms is almost constant. For 
finite wr, there appears an energy minimum so that the 
two vortices can form a bound pair; the value of A giving 
the minimum decreases with uj^. This binding originates 
from a tension Tj of the domain wall of the relative phase. 
The binding energy can be estimated as ^ T^d^n', for a ho- 
mogeneous system T^ = 8\^i\'^\'i'2\'^k/pT with the char- 
acteristic domain size A:"^ = (|*i||*2|/2wrPt)^''^ [2i|. 

The energy minimum in Fig. [Hb) is not ensured to be 
thcrmodynamically stable because the variational func- 
tion Eq. Q) does not describe the motion of the center of 
mass of a vortex molecule. We performed the simulation 
of the imaginary time propagation of the Gross-Pitaevskii 
equation derived from Eq. (Q , and found that for 51 = 
the center of mass spirals out toward the edge of the con- 
densates because of the drift instability |23. Therefore 
a rotation is necessary to stabilize the vortex molecule 
actually. Figure EJb) shows the wr dependence of the 
molecular separation 2diji for Q = 0.15. The variational 
analysis gives a good agreement with the numerically ob- 
tained value for 2dm. 

In a realistic situation, e.g., two-component BECs 
studied in the JILA group |3| , the couphng constants no 
longer satisfy ui = U2 — mi2- Then, there appear a lon- 
gitudinal magnetic field ci that aligns the spin along the 
z-axis and the spin-spin interaction C2 associated with 




tive oscillations of the lattice of vortex molecules, and 
new phenomena related with double-layer quantum Hall 
physics. 



FIG. 4: The upper figures show the profiles of the condensate 
density |\l/i|^ and |'I'2p and the lower ones the relative phase 
for Q = 0.80, OJR — 0.2 and (a) C2 = 20 (antiferromagnetic) , 
(b) C2 = —20 (ferromagnetic). The centers of mass of the 
molecules are linked by dashed lines. The bottom inset shows 
the distribution of the topological charge within [-3< x, y <3] 



Sz (antiferromagnetic for C2 > and ferromagnetic for 
C2 < [13) in Eq. Q. The structure of the vortex 
molecule for C2 — ±20 (ci—Q) is shown in Fig. O the 
detailed effects of those terms on the vortex molecule 
will be discussed elsewhere. Then, the axisymmetry of 
the topological charge and the velocity field of the vortex 
molecule is broken. For the antiferrmagnetic (ferromagn- 
tic) case, the anisotropic vorticity is distrubuted parallel 
(perpendicular) to the molecular polarization. 

As the rotation frequency is increased, a large number 
of vortices are nucleated and form a lattice. When wr 
is turned on, the vortices begin to form pairs and result 
in the lattice of vortex molecules. For ui — U2 = M12 
(ci = C2 = 0) the vortex molecules exhibits an ax- 
isymmetric feature, yielding a square lattice of meron- 
pairs (equivalently, a lattice of skyrmions under the ba- 
sis transformation of Ref. |2i||). As cjr increases fur- 
ther, the vortices of one component overlap completely 
with those of the other, forming a triangular lattice. 
If Ml 7^ M2 7^ ""12, however, the vortex molecule is 
anisotropic as shown in Fig. |31 so that the resultant 
state also exhibits a distorted lattice structure. Typical 
examples are shown in Fig. ^ where (a) and (b) corre- 
spond to the antiferromagnetic and ferromagnetic cases, 
respectively. From the inset, the direction of the ordering 
depends strongly on the distribution of the topological 
charge in each vortex molecule. 

In conclusion, we have shown that an internal coherent 
coupling between rotating two-component BECs bring 
about a new effective "vortex-molecular" field which fea- 
tures a (pseudo)spin texture with meron-pairs. The size 
of the vortex molecules can be controlled by tuning the 
external coupling field. This will open further interest- 
ing problems such as structural transitions and coUec- 
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